Discrete flat-band solitons in the Kagome lattice 
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We consider a model for a two-dimensional Kagome-lattice with defocusing nonlinearity, and 
show that families of localized discrete solitons may bifurcate from localized linear modes of the 
flat band with zero power threshold. Such fundamental nonlinear modes exist for arbitrarily strong 
nonlinearity, and correspond to unique configurations in the limit of zero inter-site coupling. We 
analyze their linear stability, and show that by choosing dynamical perturbations close to soft 
internal modes, a switching between solitons of different families may be obtained. In a window of 
small values of norm, a symmetry-broken localized state is found as the lowest-energy state. 
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INTRODUCTION 

Discrete nonlinear systems have developed into an im- 
portant area of research during the last years, where 
many theoretical and numerical predictions have been 
tested and proved experimentally [Iml. In particular, 
nonlinear optics [2] and cold atoms [3] have shown to 
provide excellent implementations to restudy old and new 
theory, mainly coming from solid state physics. The fab- 
rication of periodic and aperiodic structures, using dif- 
ferent experimental techniques, has shown an enormous 
flexibility, where simple as well as very complicated lat- 
tice topologies can be obtained in a simple and controlled 
way in different dimensions. 

There are general theories describing the bifurcation 
of families of nonlinear localized modes ("gap solitons") 
from linear band edges when a weak nonlinearity is added 
to some spatially periodic Hamiltonian system (see, e.g., 
Ref. and references therein). In generic situations, 
when the band edge is a non-degenerate local minimum 
or maximum with non-zero group-velocity dispersion, 
perturbation theory for weak nonlinearity yields a non- 
linear Schrodinger (NLS) equation for the slowly varying 
amplitude of solutions close to the band edge. In the 
most common case of an effective cubic (Kerr) nonlin- 
earity, one may then conclude from the properties of the 
corresponding NLS solitons that in one dimension (ID), 
localized gap solitons bifurcate from the linear band edge 
with zero power (norm), while in two dimensions (2D) the 
bifurcating solution appears at a non-zero power, result- 
ing in a generic excitation threshold for gap solitons in 
two (and also higher) dimensions [H. In both cases, the 
soliton envelope decays exponentially with a localization 
length that diverges in the linear limit. 

However, there are some particularly interesting classes 
of lattices where the above mentioned conditions are not 
fulfilled. Probably the most well-known example is the 
2D Kagome lattice, where one of the tight-binding bands 



(the lowest-energy one with the sign-conventions used be- 
low) is completely flat, and in addition it touches the ex- 
tremum of the second band at one point so that also the 
latter becomes degenerate (see, e.g., [6] and references 
therein). It is therefore an open question, which we here 
aim at resolving, if and how solitons may bifurcate also 
from such flat-band linear modes. 

The study of Kagome lattices has a long history, in par- 
ticular as a prototype system for geometrically frustrated 
magnetism (see, e.g., [7]). Recently, successful syntheses 
of artificial Kagome lattices have been reported in sev- 
eral contexts; e.g., nondiffracting Kagome lattice for light 
beams were obtained in [8] , photon-based litography was 
used to fabricate Kagome lattice structures in [9], in [Ifl] 
a Kagome optical lattice was realized for trapping ul- 
tracold atoms, and in [11] a metallic Kagome lattice was 
fabricated and a flat plasmonic band observed. Thus, the 
technology to experimentally observe physical phenom- 
ena arising from the presence of a flat dispersion band in 
a nonlinear lattice appears to be within reach. 

Effects of including interactions in Hubbard-type 
models on Kagome lattices have been discussed for 
fermionic ^] as well as bosonic [13] systems, and in 
both cases shown to result in an effective gap opening 
between the lowest-energy flat band and the second band 
at a certain filling factor for sufficiently strong interac- 
tion. These gap openings were related to a breaking of 
the translational symmetry of the ground states of the 
non-interacting lattices. 

Concerning localized structures in classical nonlinear 
Kagome lattices with defocusing nonlinearity, a number 
of such structures were described in Q. However, these 
authors focused on complex structures such as vortices 
and their properties in the limit of strong nonlinear- 
ity, and did not at all discuss the fundamental modes 
and their connections to the linear flat-band modes (in 
fact, Ref. does not even mention the existence of a 
flat linear band). More recent works discuss defect soli- 
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tons in Kagome optical lattices with saturable nonlinear- 
ity [TBI and localization of light in Kagome nanoribbons 
[16|; however, both these works considered exclusively 
the case of focusing nonlinearity with solitons bifurcating 
from the edge of the upper band, which is non-degenerate 
and therefore these solitons follow the standard NLS phe- 
nomenology in 2D with excitation threshold [1, 5]. 

We here study a Kagome lattice with cubic nonlinear- 
ity by considering the following DNLS-like equation: 



fore this band does not touch the fiat band at —2 



(1) 



where z corresponds to the normalized dynamical co- 
ordinate, 7 to an effective nonlinear cubic parameter, 
and Ufi represents the field amplitude at site n in a 
2D Kagome lattice [see Fig. [U^a)]. The coupling func- 
tion Vfi^rnUfn defines the linear interactions between 
Ufi and its nearest neighbors. Model ([1]) possesses 
two conserved quantities, the norm (power) defined as 
P = I'^nP, and the Hamiltonian (energy) defined 
as = - En {Ern VfiM^^K + <^fi) + {l/2)\un\^}. 
Unless otherwise stated, the defocusing nonlinear lattice 
is obtained by fixing 7 = — 1, and Vfi^rn = 1 foi" nearest 
neighbors and zero otherwise. 

Linear solutions (7 = 0) are obtained by solving 
model ([1]) with a stationary ansatz of the form Ufi{z) = 
Uj^ exjp (iXz). To obtain the linear spectrum (cf., e.g., 
[y, Illl4l3| ) , we first consider an infinite system and three 
fundamental sites belonging to a unit cell of the lattice 
[triangle in Fig. HJa)]. We construct the corresponding 
2D ^-vectors and obtain three different linear bands: 
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where f{kx^ky) = 1 + 2cos'^{kx/2) — 3 cos^(/Ca^/2) — 
cos'^{V3ky/2) + 2cos2(A:^/2)cos2(^/3A:^/2). Fig. [^b) 
shows a 3D plot of the band structure in the first Bril- 
louin zone. The upper and lower bands are "connected" 
at A = 1 by six Dirac points located at the vertices of the 
hexagon forming the Brillouin zone [f{kx,ky) = —1/8]. 
A sketch of the fundamental modes associated to the top 
of the upper band and to the bottom of the lower band 
are shown in Figs. [T]^c) and (d), for a finite-size lattice 
with rigid boundary conditions. The largest eigenvalue 
mode has a typical structure resembling the fundamental 
mode of any 2D system, where all sites oscillate in phase 
with a decaying amplitude due to the open boundary 
conditions. The smallest eigenvalue mode of the lower 
band possesses a structure reminiscent of a staggered 
mode; however, due to the particular geometry of the 
Kagome lattice it is not possible to construct a fully stag- 
gered mode, but a "frustrated" state like the one shown 
in Fig. [IJd). Note that, in contrast to the case with 
periodic boundary conditions, this eigenvalue is slightly 
larger than —2 for rigid boundary conditions, and there- 
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Figure 1: (Color Online) (a) Kagome lattice structure with 
a unit cell of three fundamental sites, including their interac- 
tions, (b) Band structure. Linear modes profiles at A = 3.91 
(c), and A = —1.98 (d), for a lattice of 205 sites with rigid 
boundary conditions. 

The third degenerated-flat band, located exactly at 
A = — 2, possesses as many states as the number of closed 
rings in the lattice (being infinite for an infinite system) 
. These states - called "6-peaks" or "ring" solutions - 
have 6 peaks with equal amplitude but alternating sign 
(phase), with a strictly zero background [see Fig. EJa)- 
inset]. These ring (hexagon) modes constitute "building 
blocks" for a Kagome lattice. Any linear combination of 
them will generate an exact linear stationary solution of 
the system. 

Therefore, a fundamental question concerning nonlin- 
ear solutions bifurcating from some specific linear com- 
bination of these modes arises. Let us first compute this 
ring stationary solution in the nonlinear regime. It is 
simple to show, that each ring mode of the flat band can 
be exactly continued into a nonlinear mode with exactly 
the same configuration^ only with a frequency shift; thus 
these solutions are "exact discrete compactons" [l7| . For 
any nonlinear ring mode, the frequency shift and power 
are related as P = 6(A + 2)/7. 

For a defocusing nonhnearity (7 < 0), all fundamental 
nonlinear solutions bifurcate at A = — 2 (P = 0) from 
some linear combination of these ring modes. Gener- 
ally in discrete nonlinear systems, a "one-peak solution" 
is often identified as a family of solutions approaching a 
single-site localized state in the strong- nonlinearity /weak 
coupling ( "anticontinuous" ) limit, and typically corre- 
sponds to the geometrically simplest fundamental nonlin- 
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ear mode bifurcating from the linear modes at the band 
borders [5]. For a Kagome lattice, a mode with these 
properties can be constructed by combining two neigh- 
boring rings, having one common central site. In the 
linear limit, this site will get twice the amplitude of the 
other ring sites. For non-zero nonlinearity, this - compact 
- solution will no longer be exact; instead, it will develop 
into a discrete soliton state with exponentially decay- 
ing tails, which continues smoothly to a single-site solu- 
tion at the anticontinuous limit (larger norms). There- 
fore, the effective size of this solution will drastically de- 
crease being, for some value of the norm, smaller than 
the ring solution. Therefore, we may expect an exchange 
of fundamental properties between the ring and the one- 
peak solutions if compared at a given norm (compare, 
e.g., similar features for saturable systems, where multi- 
ple changes on the effective size of fundamental solutions 
produce changes in their stability properties [l8l-[2Qj). 

We construct these two fundamental families (the one- 
peak and six-peaks solutions) and show the results in 
Fig. [21 The P — X diagram [Fig. Efa)] shows how these 
modes bifurcate from the flat band at P = 0, i.e., with- 
out excitation threshold. Moreover, these solutions are 
very localized for lower values of the norm being, there- 
fore, very unusual 2D discrete solitons originating in the 
particular topology of this lattice and its fundamental 
building blocks. [In fact. Fig. [2fa) is more similar to 
what is obtained for ID cubic lattices]. For larger norms, 
these two solutions strongly deviate, the norm content of 
the ring solution being much larger. 

In Fig. Efb) we show a AH — P diagram, for AH = 
Hi — Hring (^ represents any solution). For smaller val- 
ues of the norm, we see how the ring solution {AH = 0) 
corresponds to the ground state bifurcating from the 
linear band at zero norm. Then, we observe a cross- 
ing point at P ~ 0.43, between the one- and the six- 
peaks solutions, that indicates an exchange of stability 
properties. To calculate the linear stability of these so- 
lutions, a standard procedure is implemented [2l|: we 
obtain the linear eigenvalue spectrum {c<j^ = g} and 
compute the largest G = ^s/\\g\^^g]/2. A stable solu- 
tion corresponds to G = while G 7^ to an unstable 
one. The G — P diagram in Fig. [2fc) confirms that for 
small values of the norm, the six-peaks solution is sta- 
ble (ground state) while the one-peak mode is unstable. 
For 0.32 < P < 0.51, we observe an instability inversion 
regime with a region of bi-instability and, therefore, that 
none of the fundamental solutions correspond to minima 
in a Hamiltonian representation. Thus, in this regime a 
new stationary solution - corresponding to a minimum 
in-between them - should appear, connecting them in 
parameter space. For larger power, the one-peak solu- 
tion is always stable (minimum) while the ring mode is 
unstable. 

The solution corresponding to a minimum in the bi- 
unstable region is known as "Intermediate Solution" (IS) 
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Figure 2: (Color Online) (a) P vs A, (b) AH vs P, and (c) 
G YS P diagrams. The one-peak, ring, and intermediate so- 
lutions are shown with black, blue and red-dashed lines, re- 
spectively. Insets show profiles of the stationary modes for 
P = 0.43. 



and constitutes a symmetry-broken stationary solution 
appearing when the stab ility properties of fundamental 
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In the present 



solutions are exchanged 11^ 

case we find a stable IS [23|, |2^ in-between two unstable 
fundamental modes (in other contexts, the opposite is 



also possible [11, [H, S, [23| ) . In Fig. [2] (red-dashed line) 
we show the appearance of the IS connecting the two 
fundamental modes. In the stability diagram [Fig. Efc)] 
we observe how the IS is stable in the region where 
the two fundamental modes are unstable simultaneously. 
Fig. Efb) shows the emergence of the IS connecting the 
ring solution with the one-peak mode. In its existence 
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Figure 3: (Color Online) Vertical center of mass evolution of 
an unstable one-peak solution for P = 0.4655, kx = 0, ky = 
0.009. Insets: zoom of different profiles \ufi{z)\'^. Horizontal 
full (dashed) lines represent the one- (six-) peaks solutions. 



region, the IS possesses the smaller Hamiltonian value 
and constitutes an effective ground-state of the system. 

As we will exemplify by exploring the dynamics of 
an unstable one-peak solution and observing its evolu- 
tion, symmetry-broken ground states may appear spon- 
taneously. By inspecting the unstable internal mode of 
the one-peak solution [25|, we notice that it essentially 
corresponds to a deformation in the vertical direction, 
which is similar to the introduction of a phase gradient in 
the profile: ii^exp(z/c • n), with k a "kick" vector. Fig. [3] 
shows the result of applying a small kick in the verti- 
cal direction, showing a slow and smooth movement of a 
low-norm and very localized solution. We see how the ef- 
fective Hamiltonian (Peierls-Nabarro) potential is traced 
and the velocity of the center of mass changes accordingly 
(cf., e.g., [I8|,ll9|). The regions with largest velocity cor- 



respond to minima of the effective potential, which in 
this case correspond to the intermediate stationary solu- 
tions with profiles possessing a geometry in-between the 
one and the six-peaks solutions. The mobility of this so- 
lution generates some radiation that decreases the power 
of the moving profile. Although the solution starts in a 
norm region where the one-peak solution has a lower H- 
value than the six-peaks [see Fig. [2fb)], it rapidly moves 
to a region below the crossing point. Therefore, the solu- 
tion moves faster when passing through a Ring solution 
than when passing through a one-peak mode. In fact, 
for the small kick used, the solution is able to jump co- 
herently one complete site in the vertical direction and 
then it cannot continue traveling further due to radiation 
losses (bounce at z ~ 1500). Then, it oscillates passing 
through ring modes and intermediate solutions, without 
being able to overcome the one-peak barrier (horizontal 
full lines in Fig. [3]). Finally, it gets trapped oscillating 
around the symmetry-broken ground state. 

The possibility of moving a very localized solution 



across the lattice is certainly an important issue in differ- 
ent physical contexts. Typical 2D nonlinear cubic lattices 
do not allow mobility of highly localized excitations [26| , 
which are thought to be the key entities for controlling 
the propagation of information - in the form of waves - 
in periodical media. By choosing different phase gradi- 
ents, we were able to move one- and six-peaks profiles 
across the whole lattice. For example, by giving a verti- 
cal kick of ky = 0.25 to a one-peak mode of P = 0.502, 
we observe a translation of six unit cells in the verti- 
cal direction (the excitation stopped because of borders 
and radiation effects). While moving, solutions generate 
some radiation and the dynamics is not as soft as the 
one shown in Fig. [3l Nevertheless, coherent mobility of 
highly localized solutions is allowed due to the Kagome 
topology. 

In conclusion, we showed, using the Kagome lattice 
as example, how nonlinear localized modes can bifur- 
cate from a highly degenerate, dispersion-less linear band 
without excitation threshold. We identified two types of 
fundamental modes, which were shown to exchange their 
stability and therefore could be switched into each other 
through a symmetry-broken intermediate state, which 
constitutes the ground state around the exchange region. 
Since this scenario appears already for a weak (Kerr) 
nonlinearity, and involves states which are strongly lo- 
calized due to the flatness of the linear band, it could 
be of large relevance for practical applications. While 
for simple cubic regular DNLS lattices symmetry-broken 
ground states do not appear and stability properties are 
never exchanged between fundamental modes, the topol- 
ogy of the Kagome lattice allows the appearance of new 
solutions that can be crucial to improve the dynamical 
properties of these nonlinear systems. This opens new 
possibilities for controlling the propagation of waves in 
such systems. 
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Figure 4: The main figure shows the linear spectra g versus A of the one-site (blue dots) and one-ring (black dots) solutions 
in the regime of weak nonlinearity. Insets (a)-(c) show the profiles of the indicated linear modes. It is seen how the flat band 
spreads out due to the spatial symmetry breaking caused by the particular nonlinear localized fundamental mode excited, and 
how isolated localized linear eigenmodes bifurcate from this band. For the one-site solution there is one such mode which is 
unstable, while for the one-ring mode there are two stable (soft) internal modes close to the linear limit. One of these modes 
then becomes unstable in the bifurcation where the symmetry-broken intermediate solution is born (the second mode also 
becomes unstable slightly afterwards). Seen from the linear limit, one can see how a gap opens around zero, which is a classical 
counterpart to the gap openings due to interactions in quantum Hubbard models [12, [l3] • 



